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ABSTRACT 

We report the detection of strong, resolved emission from warm H2 in the Taffy galaxies and bridge. 
Relative to the continuum and faint PAH emission, the H2 emission is the strongest in the connecting 
bridge, approaching L(H2)/i(PAH8/im) — 0.1 between the two galaxies, where the purely rotational 
lines of H2 dominate the mid-infrared spectrum in a way very reminiscent of the group-wide shock in 
the interacting group Stephan's Quintet. The surface brightness in the 0-0 S(0) and S(l) H2 lines in 
the bridge is more than twice that observed at the center of the Stephan's Quintet shock. We observe 
a warm H2 mass of 4.2 x 10* Mq in the bridge, but taking into account the unobserved bridge area, 
the total warm mass is likely to be twice this value. We use excitation diagrams to characterize the 
warm molecular gas, finding an average surface mass of ^5 x 10^ Mq kpc~^ and typical excitation 
temperatures of 150-175 K. H2 emission is also seen in the galaxy disks, although there the emission is 
more consistent with normal star forming galaxies. We investigate several possible heating mechanisms 
for the bridge gas, but favor the conversion of kinetic energy from the head-on collision via turbulence 
and shocks as the main heating source. Since the cooling time for the warm H2 is short (~ 5000 yr), 
shocks must be permeating the molecular gas in bridge region in order to continue heating the H2. 
Subject headings: galaxies: individual (UGC 12914, UGC 12915) — galaxies: interactions — inter- 
galactic medium 



1. INTRODUCTION 

The Spitzer Space Telescope has led to tremendous 
advances in the study of galaxies and their interstel- 
lar media. One particularly fruitful area has been the 
study of the lower-level purely rotational emission lines 
of molecular hydrogen. An in t eresti ng example was the 
detection by lAppleton et al.l (|2006l ) of powerful emis- 
sion from purely rotational H2 lines in the group-wide 
shock of the strongly interacting group Stephan's Quin- 
tet (SQ). Spectral mapping reveals that the emission is 
distributed along the whole shock and exceeds the X-ray 
power by a factor ^ 3, indicating that H2 line emission 
is an important coo ling process in the post-shock region 
(|Cluver et alllMoh . 

Powerful H2 emission has also been detected in many 
oth er environments, including lumino us infrared galax- 



l galax.- 

ies (iLutz et al.ll2003D. radio gala xies ()Ogle et al.l I2007L 
I2010t ). A GNs (iRoussel et al.|l2007D. and cool cluste r core 
galaxies (jEgami et al.l 120061 : iDonahue et all 120 lit ). In 
another case H2 is detected in the a pparent wake of a 
bow- shock in the cluster Abell 3627 fSivana ndam et all 
|2009() . The class of strong molecular hydrog en emission 
galaxies (called MOHEGs; lOgle et al.l [20071 ) is defined 
by the strength of the H2 0-0 S(0)-S(3) lines relative to 
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that of polycyclic aromatic hydrocarbons (P AHs), with 
MOHEGs having L(H2)/L(PAH8/zm) > 0.04 (jOgle et al.l 
120101 ). The mechanism that powers the H2 emission in 
MOHEGS is not fully understood, in part because the 
emission regions are usually unresolved, but in the ra- 
dio g alaxies is likely t o be jet-driven shocks ( Ogle et aLl 
l2010tlNesvadba et al.ll2010[ ). The SQ shock is exceptional 
because it is spatially extended over ~30 kpc and asso- 
ciated with a well-studied intergalactic shock, which has 
led to a clear picture of the heating mechanism as dissi- 
pation of the kinetic energy associated with the collision 
(IGuillardet~aI1[2009l) . 

Here we report the detection of another system con- 
taining strong, resolved H2 emission in the bridge be- 
tween the interacting galaxies UGC 12914/5, also known 
as the "Taffy" galaxies because of the e xtended radio 
emiss ion stretched between the two disks (C ondon et al.l 
|1993( ). The emission closely resembles that observed in 
SQ, with the H2 lines featuring prominently in the mid- 
IR spectrum. 

The Taffy system is believed to be the result of a nearly 
head-on collision between two disk galaxies. The galax- 
ies, UGC 12914 and 12915, have heliocentric velocities 
of 4371 ± 8 and 4336 ± 7 km s~^ respectively, and we 
adopt a distance of 60 Mpc based on a Hubble constant of 
72 km s~^ Mpc~^. The disks are ~12 kpc apart and are 
believed to be separating mainly in the plane of the sky 
at a velocity of 450 km s~^ (jCondon et al.lll993f ). This es- 
timate is based on assumed masses and a parabolic orbit, 
and may be an upper limit since it does not take into ac- 
count dynamical friction resulting from the overlapping 
dark matter halos. Among the notable features of the 
system is the gas-rich bridge, which conta ins about 25% 
of the H 1 gas in the system (jCondon et al. 1993). The 
bridge is also host to a significant amount of molecular 
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gas ([Smith fc Strucld[200ll: iBraine et al.|[200l [Gao et alj 



I2OO3I) and cold dust (fzhu et al.ll2007D . A warm dust com- 
ponent is_a££arentlyjieated by UV radiation from the 
disks (iJarrett et alJ Il999[ ) . Ha images reveal that the 
bridge has little star formation outside of a single large 
H II region near UGC 12915 (Bushouse fc Werner 1990) . 

When two galaxy disks collide nearly face-on at high 
velocity, we expect widespread collisions between dif- 
fuse atomic clo uds, which have a large covering facto r 
in galaxy disks (jLerov et all 120081 : iSanchez et al.ll2010f ). 
In the collision, dense molecular clouds in one galaxy will 
ram through the more diffuse atomic gas of the second 
galaxy. The coupling of the generally smaller molecular 
clouds to larger areas via magnetic fields will increase 
the drag on the galaxies, enhancing the interaction. The 
net result of these processes is that the atomic gas will 
be splashed out into the bridge with a range of trans- 
verse velocities with values extend ing up to the rela- 
tive collision velocity ()Strucklll997[ ). Some of this gas 
will fall promptly back into both disks, and some will 
be stretched between them. Large molecular clouds will 
tend to stay close to their parent disks, which partially 
explains the lack of star formation in the bridge. The gas 
splashed into the bridge will also have a broad range of 
angular momenta, resulting in continuing cloud collisions 
that produce shocks and turbulence. 

This paper is primarily aimed at understanding the 
distribution and excitation of warm molecular hydrogen 
in the Taffy galaxies and bridge, and is organized as fol- 
lows. In Section [2] we describe the observations and data 
reduction. In Section [3] we present the mid-IR spectra, 
estimate the cooling rate, and show how the H2 emission 
is distributed in the Taffy system. We also determine the 
properties of the warm H2- In Section |4] we discuss pos- 
sible heating mechanisms for the warm H2 , and compare 
with the II2 mass inferred from previous CO observa- 
tions. We briefly summarize our findings in Section [SJ 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. IRS spectra 

We obtained spectra of the UGC 12914/5 system 
from the Spitzer public archives (Program ID: 21, PI: 
J. Houck). The observations were made wi th the In- 
frared Spectrograph (IRS; Houck et al.ll2004D o n board 
the Spitzer Space Telescope (jWerner et al.ll2004l ) on 2005 
July 8 and 10 using the Short-Low (SL; 5.2-14.5 /im), 
Long-Low (LL; 14.0-28.0 /xm), Short-High (SH; 9.9- 
19.6 /im), and Long-High (LH; 18.7-37.2 /im) modules. 
The positions of the slits for all four modules are shown 
in Figure[TJ The AORs were designed to target the nuclei 
and selected regions of the galaxies, not specifically the 
bridge. However, we were able to exploit the partially 
overlapping nature of the low resolution slits to make 
sparse maps of the region. 

The observations were initially processed by Spitzer 
Science Center (SSC) pipeline version S18.7.0. The start- 
ing point of our analysis were the Basic Calibrated Data 
(BCD) frames from the science pipeline. The BCD 
frames are a high-level calibrated data product which 
represent slope-fitted data values in electron s~^ derived 
from corrected data. Corrections include saturation, 
non-linearity, detector droop and stray-light corrections. 
BCD frames associated with each target position were 



coadded and background subtracted using dedicated off 
observations. In the SL and LL modules the dedicated 
off observations were combined with spectral orders far 
from the galaxies. Rogue pixels, especially common in 
the LH module due to cosmic-ray activation of pixels, 
were identified manually by blinking the two nod posi- 
tions against each other to help distinguish rogue pixels 
from real data. The rogue pixels were then replaced by 
suitable averaged values from adjacent pixels using cus- 
tomized software. 

Spectra were obtained at two nod positions for each 
target, with the target placed at the 1/3 and 2/3 po- 
sitions in the sht. The CUBISM software ()Smith et al.l 
l2007al) was used to construct partial spectral maps along 
the LL slits, allowing extractions from numerous posi- 
tions within the system, labeled A-U as shown in Fig- 
ure [2] The extraction regions had angular size 10'.'15 x 
10'.'15. The SL slits are oriented orthogonally and par- 
tially overlap in regions C and J. In both cases, the SL 
spectra covered about 1/3 of the area of the LL, and were 
scaled up to compensate for this difference. This rescal- 
ing assumes that the emission is uniform over the aper- 
ture, which introduces uncertainty into measurements of 
the SL lines in these regions. We therefore use these 
lines only as a check on our excitation diagrams (see 
Section [SH]). Spectra from the LL and SL slits will be 
referred to herafter as low-res spectra (resolving power 
i? between 57-127). 

The high-resolution spectra (hereafter hi-res with R = 
600) were extracted with the SSC software SPICE, us- 
ing the point source calibration for the galaxies and ex- 
tended source calibration for the bridge. The choice of 
calibrations involves some uncertainty, since the targets 
are neither fully resolved point sources nor perfectly uni- 
form extended sources. In the worst case, the difference 
between a point source and extended source extraction 
(a wavelength-dependent slit-loss correction factor) can 
lead to ^40% difference in extracted line fluxes for the 
lines observed. 

Line fluxes in each extraction regi on were measured us - 
ing the SMART software package (Higd on et aI1l2004[ ). 
In regions C and J, the partial overlap with the SL mod- 
ule expanded the wavelength coverage. In these cases 
we meas ured the line fluxes using both SMART and 
PAHFIT (|Smith et al.ll2007b[ ). PAHFIT fits known lines, 
bands, and dust contiua to a stitched version of the low- 
res spectra. Similar line fiuxes were obtained with both 
methods; results are presented in Table [1] We note a 
slight positive offset in the measured line fiuxes in region 
J by PAHFIT relative to SMART, but the differences 
are only at the 2a level. The effect is more pronounced 
in the 28 /tm line. As we point out in the footnote to 
Table [U the LLl spectrum from which this was derived 
showed some inconsistency in its continuum level com- 
pared with LL2, so this line is uncertain no matter how 
it is measured. 

The hi-res line fiuxes were measured by single Gaussian 
fitting using SMART. The LH slit covers an area 22'.'3 x 
ll'.'l. The SH slit covers a smaller area of 11'.'3 x 4'.'7. 
The SH spectra were scaled at each nod position so that 
the continuum matched between the two modules. The 
scaling factors were similar for both nods. The line fluxes 
and mean of the scaling factors for the two nods are 
presented in Table [1] 
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2.2. IRAC and MIPS images 

Images of tfie system obtained w itli tlie Infrared Array 
Camera fIRAC: lFazio et al.ll20M ) an d Multiband Ima, 
ing Pliotometer for Spitzer (MIPS; iRieke et all [20a 
were available from the Spitzer public archives (Program 
ID: 21, PI: J. Houck). These data were processed by SSC 
pipeline version S18.7.0, and were of sufficient quality 
that further processing was not required. 

Photometric measurements were made in each of the 
LL bridge extraction regions using square apertures with 
the IRAF task polyphot. The CUBISM spectral extrac- 
tions use extended source calibration, so the photometric 
measurements were not aperture-corrected. The IRAC 
3.6, 4.5, 5.8, and 8.0 /xm images were Gaussian-convolved 
to match the resolution of the MIPS 24 /im image, which 
is comparable to the resolution at the II2 S(0) 28.22 /xm 
line. For comparison with the hi-res spectra, we also 
measured the images using a set of boxes matched to 
each of the LH nod positions. 

3. RESULTS 
3.1. Spectra 

The hi-res spectra of the nuclei and the bridge are 
shown in Figure [3] The strength of the H2 lines is very 
striking, especially in the bridge, where the 0-0 S(0) and 
S(l) lines are stronger than all other lines except the 
[Si II]34.82 /im line. The S(l) line is similarly powerful 
in the UGC 12914 nucleus. The II2 line fluxes are pre- 
sented at the top of Table [I] while the fine structure line 
fluxes are shown in Table [5J We will further quantify the 
strength of the II2 emission in Section 14.11 

The low-res spectra show that the emission in the 
bridge is not confined to the position of the high reso- 
lution bridge/arm slit which, as shown in Figure [U also 
includes the UGC 12914 spiral arm. The H2 S(l) line 
is measured in regions A-S. The two regions where it 
could not be measured (T and U) are near the southern 
end of UGC 12914 and lie outside the app arent extent of 
the br idge, based on the radio contours of iCondon et al.l 
(jl993[ ). Spectra from all of the regions are shown in 
Figures |4] and [5] The red diamonds show the contin- 
uum flux densities measured from IRAC and MIPS for 
the same extraction regions. Note that in almost all 
cases these broadband measurements show good agree- 
ment with the spectra, indicating that the relative scaling 
made between spectral orders was appropriate. In only 
one case (region O) do the IRAC fluxes disagree with 
the extracted spectra. However, since no lines were de- 
tected in that position, the discrepancy does not effect 
our analysis of the H2 . 

3.2. H2 distribution 

We now present two alternative ways of demonstrating 
that the H2 emission, as measured in the low-res spec- 
tra, extends between the galaxies in the bridge. Figure IH 
shows cross sections of the LLl and LL2 slits which cut 
through the system along three slices, which we call posi- 
tions 1-3. The locations of the galaxies at each position 
are shown in the top panels. No significant emission was 
detected at position 4, so it is not shown. 

Several important results emerge from the cross- 
sectional analysis. First, both the [S III]33.48/im and the 
continuum emission (measured at wavelengths of 30 /Ltni 



and 17.4 /im) are concentrated in the galaxies, with little 
originating in the bridge. This contrasts with the H2 S(0) 
and S(l) lines, which show significant emission between 
the galaxies, especially at position 2. The H2 S(l) line 
is even stronger in the bridge compared to the galaxies, 
and actually drops significantly just south of UGC 12915. 
The right panels of Figure [HI show that the emission is 
much weaker outside of UGC 12914 at position 3 than 
at position 2 (see also Figures |4] and [5]) , indicating that 
either the gas density or temperature drops significantly 
south of a line connecting the nuclei. 

The extended nature of the H2 emission in the bridge 
can also be demonstrated using sparse spectral maps. 
Although the Taffy system was not fully mapped, suf- 
ficient long-slit low-res spectra were obtained to allow 
the construction of sparse LL spectral maps using CU- 
BISM, which are shown in Figure [71 To provide the 
reader with a better visualization of the results, the S(0) 
and S(l) data (color image) are overlaid with 20 cm 
radio continuum (contours) from ICondon et al.l ()1993[ ) . 
We also show these radio contours overlaid on an IRAC 
8.0 /im image, and sparse maps of the [S III]18.71/im and 
[Si II]34.82/im lines. These latter lines are selected be- 
cause they have wavelengths comparable to the H2 S(l) 
and S(0) lines, and so provide a sense of the contami- 
nation level in the bridge as a result of emission from 
the galaxies. Figure [7] shows that extended H2 emission 
is present i n the bridge region , similar to the distribu- 
tion of dust (iJarrett etalll 19991 IZhu et al.ll2007[ ) and CO 
(jGao et all [20031), thou gh differing from the H I, which 
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peaks in the bridge (jCondon et al, 

3.3. Total H2 cooling rate 

The total H2 luminosity (or equivalently, the II2 line 
cooling rate) in various regions of the bridge, and an es- 
timate for the whole bridge, can be calculated from the 
II2 line fluxes. Since the spatial coverage is more com- 
plete for the low-res spectra, we restrict our discussion 
to those data, which cover the S(0) and S(l) lines only, 
with the exception of regions C and J. 

To determine the total 0-0 S(0) and S(l) line lumi- 
nosity in the bridge, we exclude the regions obviously on 
the disk of the galaxies (A, G, H, P, and T), and region 
F (which contains the star forming knot), as well as re- 
gions very near the galaxies, where contamination from 
the disk is likely to be largest (B and I). This leaves re- 
gions C, D, E, J, K, L, M, Q, and R as bridge spectra. 
The total flux from these regions is 5.8 x 10~^^ W m~^, 
giving a total luminosity of 1.16x10^^ W over a pro- 
jected area of 78.5 kpc^. The mean surface brightness is 
~1.5 X 10^2 W kpc-2[] 

To estimate the total H2 luminosity, we must also es- 
timate the spatial extent of the emission. We assume 
that the II2 emission is roughly bounded by the extent 

^ For the hi-res region that lies in the bridge (called 
"bridge/arm" in Table 1) we note that the surface brightness of 
the H2 0-0 S(0) and S(l) lines, averaged over the two LH nod po- 
sitions (covering an area of 29'.'7 X ll'.'l) is 1.0 X 10^^ W kpc~^. 
This is lower than the surface brightness seen in region C mea- 
sured in the low-res spectra. Although we do not have a complete 
explanation for the discrepancy, the result may indicate that the 
emission is clumpy on the scale of the LH slit. We note that surface 
brightness variations comparable to this difference are seen on the 
scale of the low- res apertures (~10"), which may partly explain 
the difference. 
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of the radio emitting bridge. We define this area as being 
bounded by the extent of the radio continuu m (down to 
a level of 0.1 mJy in the lCondon et allll993l 20 cm emis- 
sion map) assuming reasonable inner edges to the galaxy 
disks. This provides a bridge area of 170 kpc^, and gives 
a total bridge Ha luminosity L(H2) = 2.6(±0.7) x lO^^ W, 
assuming the properties of the unobserved bridge are 
similar to those of the observed component. In regions 
C and J where the SL data allows the addition of the 
S(2) and S(3) lines, the extra contributions amount to 
an increase of 28% in the luminosity. Our estimate of 
the total warm H2 luminosity is thus likely to be a lower 
limit. 

3.4. Excitation diagram and H2 mass surface densities 

Although we do not have full spectral coverage of the 
whole bridge, we can explore the variation of H2 prop- 
erties as a function of position in the bridge where we 
have useful data. The extraction regions A-U, defined in 
Figure[2l are generally restricted to the LL module of the 
IRS, and thus cover the only 0-0 S(0) and S(l) lines. As 
previously discussed, regions C and J also provide infor- 
mation for the S(2) and S(3) lines at 12.3 and 9.66 /im. 
These measurements are sufficient to allow us to make 
a preliminary exploration of the excitation properties of 
the warm H2 across the bridge, subject to several con- 
straints and assumptions discussed below. 

We have constructed II2 excitation diagrams for the re- 
gions A-U. These diagrams plot the column density Nu 
of H2 in the upper level of each transition, normalized 
by its statistical weight, versu s the upper level energy E^ 
(e.g., iRigopoulou et al.l 120021 ). which were derived from 
the measured fluxes assuming local thermodynamic equi- 
librium (LTE) for each position. Since most of the exci- 
tation diagrams consist of only two points, corresponding 
to the 0-0 S(0) and S(l) lines, we do not show them here, 
but rather present in Table [3] the results for a single- 
temperature fit through the two points in each region. 
These fits provide a baseline measurement, and allow us 
to explore how the ratio of the S(0) and S(l) lines might 
vary as a function of position in the system. For regions 
N, O, and S only an upper limit was obtained for the 
S(0) line, so no temperature was derived. No II2 lines 
were detected in region U. 

The results of this analysis are presented in Table [31 
Most of the regions in the bridge have temperatures 
(represented by the slope in the excitation diagram) in 
the range 157-175 K, with only two regions (Q and R) 
having lower temperatures of 130-133 K. The nuclei of 
both galaxies (regions G and H) show a higher tem- 
perature. The UGC 12914 nucleus has a H2 tempera- 
ture of 195±12 K, some 30 K warmer than the average 
bridge region. The warm H2 surface density derived from 
these measurements (also given in Table [3]) range from 
lO"^ Mq kpc~^ at position F, which contains a star form- 
ing knot, to the lowest value of 2.4 x 10^ Mq kpc^^ at 
region M. The mid-bridge regions D, K, and R have sim- 
ilar H2 surface densities of 5.0-6.6x10^ Mq kpc~^. We 
again select regions C, D, E, J, K, L, M, Q, and R as 
representing the bridge. The mean and median value of 
the surface density over regions between the galaxies is 
5.3 ±0.5 and 5.5 x 10^ Mq kpc~^, respectively. Based on 
the cooling rate from the same regions (see Section [ 



and an assumed average temperature of 166 K, the av- 
erage thermal cooling time for this warm mass density is 
~5000 yr, which is very short compared with the bridge 
formation time. 

Summing the observed warm H2 masses over these 
representative bridge regions (78.5 kpc^ in area) yields 
an observed mass of 4.2 x 10^ Mq. We estimate that 
the total warm H2 mass for the whole bridge area of 
~ 170 kpc^ could be as large as 9(-h2/-5)xl0* Mq if 
the entire bridge has similar properties to the observed 
bridge regions. The lower limit comes from the extreme 
assumption that no new H2 is detected in the unobserved 
area. 

We have made several simplifying assumptions in es- 
timating the H2 surface density from these data. First, 
we have assumed a thermalized equilibrium value for the 
ortho-to-para ratio. For temperatures above 300 K, this 
ratio is 3, but for temperatures in the range 130-180 K 
this value varies froni 2.3-2 .8 respectively (see Equation 
4 of IWilgenbus et ahl 120001) . Under most circumstances 
we might expect the lower- J transitions to be in equilib- 
rium. Normally, the ortho-to-para ratio would be inves- 
tigated by looking for systematic differences in the odd 
and even transitions in the excitation diagram. However, 
because we have so few points (in most cases only two) we 
must simply assume equilibrium values. This obviously 
introduces an unknown uncertainty in the final derived 
properties of the gas. Deviations from LTE would lead 
to uncertainties in the assumed excitation temperature 
of the gas and thus the final total column densities. In 
regions C and J where four lines are detected (discussed 
below) we do not see any obvious deviations in the odd 
and even values for N/g versus upper-level energy, but 
this trend may not hold in all bridge regions. 

A second assumption is that the 0-0 S(0) and S(l) 
transitions can be fitted by a single temperature com- 
ponent. Excitation diagrams for H2 observations are 
usually fitted by se veral temperature components (e.g., 
iRoussel et al.ll200"7l ). To explore this assumption we use 
the low-res regions C and J, where we were able to mea- 
sure the S(2) and S(3) lines, to investigate how this as- 
sumption affects the results. The fits for these regions 
are shown in Figure [8l while the gas parameters are pre- 
sented in Table [3] We find that in addition to a cool com- 
ponent, a warmer component (^430-440 K) is needed 
to explain the S(2) and S(3) measurements. However, 
the effect on the temperature of the coolest component, 
which dominates the mass surface density, is quite small. 

A two-temperature fit can also be derived for hi-res 
bridge/arm region, since the hi-res data includes the S(2) 
line. Wc find Ti — 102 K for the coolest component, and 
T2 = 310 K for the warmer one (Figure [HI right panel). 
However, with only three observed points there is signif- 
icant degeneracy between the temperature and the col- 
umn density, so the results should be taken as very ap- 
proximate and poorly constrained. Adopting the above 
temperatures, the column density of 1.5 x 10^ AIq kpc~^ 
in the cooler component would be 2.5 times the value de- 
rived from the low-res data, and is probably unrealistic. 
We consider them as strict upper limits only. 

Table [3] also tabulates the H2 line fluxes for the 0- 
S(0)-S(2) lines for the nucleus and outer disk positions 
of UGC 12914, and the UGC 12915 nucleus as measured 
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from the hi-res spectra. Figure [S] shows the fits to the 
excitation diagrams for the positions on the galaxies in 
the left panel. 

4. DISCUSSION 

4.1. H2 heating sources 

4.1.1. PDR heating 

To quantify the strength of the H2 emission, we com- 
pare the total flux in the H2 0-0 S(0) and S(l) lines 
with the PAH emission in the IRAC 8.0 /im band. The 
PAH8/zm strength is estimated using the method of 
iHelou et al.l (|2004[ ). in which the flux from the starlight- 
dominated 3.6 /xm band, scaled by a factor 0.232, is 
subtra cted from the flu x in the 8.0 ^m band. Fol- 
lowing lOgle et all (|2010t ). we plot the luminosity ratio 
L(H2)/L(PAH8/im) against the 24 /im luminosity z^L^, 
for both the low- and hi-res extraction regions in Fig- 
ure 1101 The 24 /xm luminosity is determined from the 
MIPS measurements for all of the data points, except 
for the hi-res apertures. For the hi-res apertures associ- 
ated with the galaxies, the MIPS measurements are very 
sensitive to the exact centering of the extraction box, 
and so we self-consistently used the average flux over 
the range 22.5-25.0 /im from the spectrum to form the 
24 /xm continuum. This also ensures that any bias in- 
troduced in the extractions by assuming the lines come 
from a point source rather than an extended source do 
not somehow influence the results since the slit correc- 
tion factor for point sources would disappear in the ratio. 
In addition to the Taffy regions, we show star forming 
galaxies, LINERS, a nd Seyferts from the SINGS sample 
(iRoussel et al.ll2007D . We also show the SQ shock sub- 
region CGluve r et al.l 120101 ) and Arp 143 knot G, which 
had the highest L(H2 )/£(PAH8/im) rati o of any part of 
the Arp 143 system (JBeirao et alll2009| ). F{PAR8^m) 
has been determined in the same way for all data points. 

In Figure [TO] we have distinguished between spectra 
taken on the bright parts of the galaxy disks and spec- 
tra taken between the galaxies in the bridge. The spec- 
tra centered close to bright disk regions (or the nu- 
clei themselves) show L(H2)/i(PAH8/xm) ratios typical 
of star-forming galaxies and likely associated with pho- 
todissociation regions (PDRs) around young stars. The 
L(H2)/i(PAH8/im) ratio in the Taffy bridge regions are 
much larger than in the bright disks, and are generally 
higher than those of the SINGS AGNs. We note that 
Region F, which includes the giant H II region, shows a 
higher L(H2)/i(PAH8/xm) ratio than is typical of star 
forming galaxies. This is not surprising, since it lies 
in the bridge and may have both bridge and PDR-type 
heating present along the line of sight. 

There is generally good agreement between the high 
and low resolution data in the emission regions within the 
galaxies. The hi-res spectra of the UGC 12915 nucleus 
cover the same part of the system as region G, the UGC 
12914 nucleus is comparable to region H, the northern 
clump of UGC 12914 is comparable to region A, and the 
UGC 12914 southern knot is not far from region P. These 
regions have L(Il2)/i(PAH8/im) ratios on the high end 
of the star forming galaxy distribution which, along with 
their positions within the galaxies, suggests that the ex- 
citation mechanism may be the same as that of the star 
forming (not AGN dominated) SINGS galaxies, namely 



PDR heating. 

A similar conclusion is reached by exa mining Fig- 
ure[TTl which shows the ratio L(H2)/L24, as in I Ogle et al.l 
([2010). In addition to the SQ shock s ub-region and Ar p 
143 knot G, we also plot here 3C 326 N (|Ogleet al.l2 007y 
Emission at 24 /xm is mostly due to warm dust heated 
by young stars. The relative H2 emission in the Taffy 
bridge regions is considerably stronger than that in the 
star forming galaxies, and is also stronger than most of 
the AGNs. The bridge regions are generally comparable 
to 3C 326 N and Arp 143 knot G, but somewhat weaker 
than the SQ shock sub-region. 

The possibility of PDR heating can also be inves- 
tigated by coinparin g to the Meudon PDR models of 
iLe Petit et al.l ()2006[ ). The models use the flux ratio of 
the H2 S(0)-S(1) lines to the CO(l-O) line to measure 
the heating rate in the in olecular gas (jNesvadba et al.l 
20101 iGuillard et alll2012| ). We use the parameters of 



Habart et al.l (120111 ) to produce four different models 
with riH = 100 and 1000 cm~'^ and UV radiation scal- 
ing factor Guv = 1 a-nd 10. The radiation scaling: 
factor is deflned relative to the UV field of Habin^ 
(11968) . The models predict i^(H2)/F(C0) w 20-60 
for PDR heating. For co mparison w i th th e models, we 
use the CO(l-O) data of iGao et all (|2003[ ). which give 
i^(CO) = 1.2 X 10^1^ W m-2 over the entire bridge. Us- 
ing the total H2 bridge flux of 5.8 x 10~^^ W m^^, we 
find F(Il2)/i^(C0) = 460, about an order of magnitude 
above the predicted levels for PDR excitation. 

The two galaxies are not exceptionally strong H2 emit- 
ters compared to the SINGS star forming galaxies, with 
L(H2)/i(PAH8/im) ratios only slightly above the most 
Il2-bright star forming galaxies. In contrast to most of 
the bridge, the two galaxies both host star formation. 
UGC 12915 seems to be in an early starburst phase 
parrett et al.| [l"999). so much of the H2 emission is likely 
to have a PDR origin. Star formation i n UGC 12914 
does not appear to be as strongly enhanced (jJarrett et al.l 
Il999f ) , but the II2 emission is consistent with star forming 
galaxies. 

For comparison with the bridge, we compare the 
i^(H2)/F(C0) ratios of UGC 12915 nucleus (region G) 
and UGC 129 14 South (^ r egion P). Using the CO(l-O) 
line profiles of IGao et al.l (|2003f ) and rescaling to match 
the areas of bridge region apertures G and P, we find 
CO fluxes i^(CO) = 9.4 x IQ-i^ and 3.3 x IQ-i^ W m'^. 
This gives F(H2)/i^(C0) = 78 and 98 in UGC 12915 
and UGC 12914 South, respectively. Both galaxies have 
stronger II2 emission than expected for PDR heating, but 
are within a factor of two of the models. It seems plausi- 
ble that PDR heating is a signiflcant source of heating in 
the galaxies, but there are probably additional heating 
sources as well. More significantly, the line ratios in the 
galaxies are considerably smaller than that in the bridge, 
highlighting the likely difference in the dominant heating 
mechanism. 

4.1.2. Cosmic ray heating 

Given the relatively strong radio emission from the 
bridge region, we consider whether cosmic ray (CR) heat- 
ing could be responsible for the excitation of the warm 
II2 emission. This process involves low-energy cosmic 
rays (1-10 MeV), which ionize some of the gas. The 
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resulting primary and secondary electro ns heat the gas, 
which excites the H2 through collisions (|Dalgarno et alj 
I1999D . Radio continuum observations of the bridge pro- 
vide some clues as to the energy density of cosmic rays. 
We caution however that the CRs that are expected to 
excite H2 are at the very low energy tail of the CRs that 
give rise to the radio continuum emission. 

The VLA 20 cm radio conti nuum flux in the mid-bridge 
is ^ 1.1 mJy per 5" beam (jCondon et al.| [r993). We 
assume a power law spectral distribution of the form 
i/~" with a = 1.2 over the range 1.49-4.96 GHz, and a 
plasma depth of 12 kpc (roughly the scale of the bridge). 
The minimum equipartition magnetic field is found to 
be 6.4 yuG fe.g.. lGovoni fc Ferettill200^ for a lower spec- 
tral cut-of f at 10 MHz. This is close to the value ob- 
tained by ICondon et al.l ([l993) under similar assump- 
tions. The corresponding magnetic energy density is 
3.8 X 10"^'^ J m~'^. This energy density is uncertain to 
within at least a factor of two because of uncertainties in 
the value of the proton/electron energy ratio (assumed 
unity), as well as uncertainties in the volume filling fac- 
tor of the magnetic field (also assumed to be unity) and 
the depth of the plasma column. 

If the energy density in the cosmic rays < C/cR > ~ 
< Ub >min, then we can estimate the power available 
for heating Phcat — < Uqr > x ?7/t, where i] is the 
efficiency of the CR heating of H2, and r is the char- 
acteristic deposition timescale. For the central bridge, 
assuming a depth of 12 kpc, this corresponds to a lumi- 
nosity surface density Lqr — 4.5 x 10^^ [v/ti) W kpc^^, 
where T7 is the deposition timescale in units of 10 Myr, 
the approximate expansion timescale of the bridge. We 
estimate 77 by approximating the energy deposited by a 
typical MeV cosmic ray from its specific stopping power 
measured in MeV per unit column density. A value of 
3.5 MeV g~ ^ cm~^ is quoted for typic al MeV cosmic rays 
in the ISM (jYusef-Zadeh et al.|[2007l ). where the column 
density is that of the target material, in this case a mix- 
ture of molecular hydrogen and H I. With H2 column 
densities in the Taffy bridge of ^2 x 10^° molecule cm~^, 
the stopping power liberated by the passage of these CRs 
through this medium would be of order 3.5 x 10^^ MeV 
per particle. Even allowing for the existence of more cold 
H2 material, this implies 77 < 0.01. In order to power the 
H2 for the lifetime of the bridge {tj — 1) with a mea- 
sured H2 surface brightness of ~1.5 x 10"^^ W kpc~^, a 
cosmic ray energy density of at least 100 times that de- 
rived from equipartition arguments would be required. 
Alternatively, CRs could heat the H2 if we are observing 
the system during a deposition burst (ry ~ 0.01). Such a 
process would have to be global, perhaps taking the form 
of a sudden injection of CRs streaming from the galactic 
disks, but this seems unlikely. 

Another approach to evaluating the infiuence of cos- 
mic rays, which is independent of assumptions about 
equipartition, is to estimate the ionization rate CcR 
needed to balance the H2 line cooling if CRs were 
the primary heating source. The adopted H2 surface 
brightness of 1.5x10"^^ W kpc~^ gives a cooling rate 
'--^ 5 X 10"'^^ W molecule"^, assuming an average mass 
density of 5.3 x 10^ Mq kpc~^ over the bridge. The heat- 
ing rate by CR ionization has bee n estimated by various 
authors. lYusef-Zadeh et all (J2007D estimate the CR heat- 



ing to be 2 X (4 X IO^^^Chs) W molecule"^, where Chs 
is the H2 ionization rate. Balancing H2 cooling with CR 
heating under these assumptions would require an ion- 
ization rate of ^1 0~^^ s~^. Under the s a me as sumptions, 
IQgle et al.l ()2010l ) and 'Nesv adba et al.l ([2010) infer sim- 
ilar values for the ionization rate in the MOHEG radio 
galaxy 3C 326, where shocks were implicated. 

This ionization rate is significantly higher than that 
measured in the molecular clouds in the Galactic Center 
(jOka et al.ll2005D . Thus, on purely comparative grounds, 
a CR ionization rate in the Taffy bridge would again 
require an unusually high CR energy density, perhaps 10 
times that of the unusual Galactic Center regions and 
~100 times the galactic neighborhood. Based on these 
arguments we can conclude that cosmic rays are unlikely 
to be responsible for heating the H2 in the Taffy bridge. 

4.1.3. Heating by magnetic reconnection 

This process is closely related to cosmic ray heating. 
In the previous subsection we estimated that the mag- 
netic energy density in the bridge plasma was about 
3.8 X 10"^'^ J m~'^. Again assuming a bridge depth of 
about 12 kpc, the energy column density is ~ 1.4 x 
10^^ J kpc^^. Supposing that this magnetic energy is ex- 
tracted on the bridge expansion timescale of about 10^ yr 
(jCondon et al.|[l993l ). we expect a surface luminosity of 
about. Free = 4:A x lO^Vy W kpc'^. 

This is comparable to the corresponding cosmic ray 
surface luminosity, as would be expected from the 
equipartition approximation of the previous subsection. 
Like the cosmic ray luminosity, this estimate must be 
corrected for an efficiency factor, representing the frac- 
tion of the reconnection energy used to excite H2 . This 
should include both direct excitation from the radiation 
from reconnection regions, and indirect processes, like 
broader ambient heating from reconnection. In either 
case we expect a very low efficiency, since the plasma will 
be transparent to most wavebands of the broad radiation 
spectrum produced directly or indirectly in reconnection 
events. 

Note that the surface luminosity estimate above as- 
sumes the extraction of all of the magnetic energy on the 
adopted timescale, and so is an overestimate. On the 
other hand, it does not account for additional field gener- 
ation in turbulent dynamos. This process may be locally 
important, but it is hard to see how it could have a global 
effect in the short time available. With the magnetic 
field estimate a bove, and mass d ensity 10^ Mq kpc^^ 
estimated from IGao et all ()2003D . the Alfven velocity 
va ~ 20 km s~^, an order of magnitude les s than the 
bridge expansion velocity of 450 km s^^ (Co ndon et al.l 
|1993[) , so even the propagation of magnetic effects must 
be very localized. To create significant power over the 
whole bridge would require some organized triggering of 
many localized reconnection events. One possible mecha- 
nism is large-scale turbulence, to be discussed in the next 
section. Turbulence could, in principal, create locally 
tangled magnetic fields which could then be stretched by 
the galaxies as they move apart, creating the conditions 
for possible reconnection. If the turbulence continued to 
be present over a long period of time, new events could 
continue to be created which might heat the H2. Thus, 
although unlikely, we cannot completely rule out mag- 
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netic reconnection as another source of H2 heating. 

Specifically, if the typical reconnection scale relative to 
the bridge size is much less than the ratio of the Alfven 
speed to the bridge expansion speed, then the reconnec- 
tion timescale is much less than the age of the bridge. 
In that case, with a short timescale for reconnection, the 
reconnection luminosity is likely to decay on a similar 
timescale, yielding less luminosity than at earlier times. 

4.1.4. Turbulence and/or shock heating 

A more likely source for H2 heating is turbulence, cre- 
ated in the wake of the collision between the two galaxies. 
We estimate turbulent heating rate Fturb using 

_ 3 2toh(w/2.36)2 

t turb — ;^ X - , 

where v is the velocity width over a length scale I, 
the factor 2.36 converts the line width to the rms veloc- 
ity of the gas, and idis = //(z;/2.3 6) is the energy dis- 
sipation timescale (e.g., iMac Lowl [l9 99: Ticlcns 2005). 
From the CO observations of iGao et al.. (,2003,), wc take 
V « 200 km s^^ over a beam of 14", or Z Ri 4 kpc. 
This yields fdis « 5 x 10^ yr. The heating rate is then 
Tturb ~ 2.5 X 10-32 W molecule-i, about 50% of the 
observed cooling rate of 5 x lO-'^^ w molecule"^. Tur- 
bulence on this scale could contribute significantly to the 
gas heating, provided that it is relatively efficient. 

A related possibility is heating by shocks. As noted 
earlier, the H2 line cooling time is very short, and in-situ 
heating by shocks or turbulence is consistent with the 
need to maintain the observed gas at ~160-170 K for a 
reasonable fraction of the bridge lifetime. 

We may obtain a crude estimate of the energy avail- 
able in shocks by examining the bulk mechanical en- 
ergy associated with the post-shock gas. The CO ve- 
locity dispersion, which should roughly track the ve- 
locity dispersion o f the cold H2, is (200/2.36) km s"^ 
(|Gao et al.ll200l . iZhu et all (pOOl found the H2 mass 
to be Mh2 ~ 1-3 X 10^ Mq, so over the 20 Myr since the 
collision the average heating rate is 1.5 x 10'^^ W. This is 
comparable to the bridge H2 luminosity of 2.5 x lO'^'* W. 

This estimate can be checked against an upper limit 
to the available power estimated using the bulk kinetic 
energy of the H I and H2 gas in the bridge . A re-analysis 
of the radio data of C ondon et al.l ()1993[ ) by IGao et ahl 
(|2003) determined that the H I mass in the bridge is 
Mhi ~ 6 X 10^ Mq. We take this mass as the pre- 
shock H I mass, and a collision velocity of 450 km s~^ 
(jCondon et al.|[l993D in the shock frame. The time since 
the collision is ^ 20 Myr, so the total available power 
from the H I is 2 x 10"^^ W, about 2 orders of magnitude 
higher than the H2 luminosity. 

Of course, not all of the available power from any heat- 
ing mechanism will be dissipated via H2 line emission. 
We note that H2 will not be the only line-coolant in the 
shock, although it could be a significant one. For ex- 
ample, in our modeling of the similar shock-heating of 
the Stephan's Quintet system, we estimate that 75% of 
the line cooling in the molecular shocks is from H2, with 
the majority of the other cooling through H2O, CO and 
[O I]/xm line emission. The details of the cooling channel 
wou ld depend s omewhat on the mix of C- and J-shocks 
(see iFlower (SFPincau dcs Forcts,i2010.) . 



Other, higher velocity shocks which might be present 
could also dissipate significant energy through mecha- 
nisms such as Ly-a emission and other UV coolants, as 
well as diffuse X-ray gas heated in the collision. There 
is not currently any X-ray data available for the Taffy 
that could be used to determine the significance of X- 
ray cooling in the system. In Stephan's Quintet, the to- 
tal X-ray lumin osity is about 3 tim es lower than the H2 
line luminosity (jCluver et al.ll2010l) . but it is not clear 
that the same fraction should apply to the Taffy. The 
presence of shock s is required by the recent model of 
iLisenfeld fc Volkl (|2010( ). In the model, strong shocks in 
the Taffy bridge accelerate charged particles to ultrarel- 
ativistic energies, which then produce the synchrotron 
emission observed in the radio continuum. 

We conclude that the largest heating s ource is proba- 
bly shocks, as in the model presented bv iGuillard et all 
(120091 which was developed to explain the powerful H2 
emission detected in the giant intergalactic shock wave 
in Stephan's Quintet. The essential ingredients of this 
model are a high-speed shock wave that is driven into a 
multi-phase medium, leading to the collapse and shock- 
heating of denser material. Molecules form rapidly on 
dust grains which survive in the denser clumps of mate- 
rial, and strong H2 emission results. Models fitted to the 
II2 excitation diagrams in Stephan's Quintet show that 
the emission is consistent with several low and medium 
velocity C-shocks, although more destructive J-shocks 
cannot be ruled out. Our observations do not provide 
enough points on the excitation diagram to justify exten- 
sive shock-modeling, but the similarities with the spectra 
seen in the Taffy bridge (powerful H2 with only low- 
excitation metal lines) and the approximately similar 
low- J temperature fits suggest a common origin for the 
II2 excitation. 

By analogy with Stephan's Quintet, a distribution of 
shock velocities, related to the dens ity distribution in th e 
shocked gas, is likely to be present (jGuillard et al.|[2009l) . 
In the gas that is less dense than the molecular gas, shock 
velocities > 100 km a~^ are likely, which will ionize the 
medium. These shock velocities can be estimated us- 
ing the observed ratio [Ne III]15.56 /im/[Ne II]12.81 /im, 
which we obtain from the hi-res bridge/arm region where 
both of these lines were observed with the SH module0 
These lines can also originate with star formation, and 
the observed line ratio [Ne III]/[Ne II] w 0.62 in the 
bridge/arm regi on is within the range typical of star 
forming regions (jDale et al.ll2006[ ). but Ha maps show 
no evid ence of significant star fo rmation in the Taffy 
bridge (jBushouse fc Wernei1ll990( ). though star forma- 
tion is likely to be present in the UGC 12914 spir al arm. 
We u se the MAPPINGS shock models of .Allen et al.l 
(|2008D for a magnetic field of 5 fiG and a pre-shock 
gas density between 0.1-100 cm"'^, as shown in Fig- 
ure [T^ The shock velocity for the Taffy bridge falls in 
the range 1 00-300 km s"'^, comparable to the SQ shock 
(jCluver et a l. 2010). 

Figures [10] and [11] may provide some insight into the 

* The [Ne II] 12.81 ^m was observed in the low-res data only 
in regions C and J. Since it falls on a separate module from the 
[Ne III]15.56 ^m line, the ratio will depend on the scaling factor 
chosen between the two modules. We therefore opt instead to 
use the hi-res bridge/arm region, where both lines fall on the SH 
module and the ratio is independent of scaling. 
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location of the strongest shocks. Figure [TUl shows a spa- 
tial trend in L(H2)/L(PAH), with values close to those 
of the galaxies at the bridge/galaxy interface growing to 
become comparable to the SQ shock near the center of 
the bridge. A similar trend is found in L{}l2)/L24, shown 
in Figure [TTJ This suggests that the shocks in the Taffy 
bridge are strongest relative to other excitation mecha- 
nisms near its center. 

4.2. Additional comparison with other systems 

It is worth making some additional comparison be- 
tween the properties of the Taffy bridge and those in the 
large-sc ale shocked filaments in the Stephan's Q uintet 
system (jAppleton et al.l l2006t iCluver et all l2010( ) . The 
mean surface brightness of the low resolution S(0) and 
S(l) lines is '--'1.5x10^^ W kpc~^, more than a factor 
of 2 larger than the equivalent emission averaged over 
the inner region of the main shock in Stephan's Quintet, 
which we estimate to be 7x10^^ W kpc^^ for the S(0) 
and S(l) lines. Thus the surface brightness in the Taffy 
is greater, reflecting a larger surface mass density of ex- 
cited gas. We can probably rule out higher excitation of 
the warm gas as the explanation because, as we show in 
Section l4?3l the ratio of warm to total H2 gas in the Taffy 
is probably lower (within the considerable uncertainties) 
than that in SQ. Thus the higher surface luminosity is 
almost certainly a result of the different quantities of gas 
being swept up in the shock in the Taffy compared with 
SQ, because of the different nature of the collisions in 
the two systems. 

Alternatively, the higher surface density of warm gas 
may be purely geometrical - the Taffy bridge could be 
deeper in the plane of the sky than the SQ shock, a re- 
sult of very different collision geometries in the two cases. 
The total gas mass in the Taffy bridge would depend on 
many complicated factors resulting from the geometry 
of the original collision, including disk impact parame- 
ters and inclination and the rotational kinematics of the 
disks. All of these factors would govern how much mass 
is splash ed into the b ridge and its projected surface den- 
sity fsee Struckl ll997l ). 

Local Luminous Infrared Galaxies (LIRGs) and Ultra- 
luminous Infrared Galaxies (ULIRGs) comprise a high 
fraction of collisional or merging systems. The 
Taffy system has a far-IR luminosity as measured by 
iSanders et all (pOOS) Lir = 10^°-9 Lq, just short of the 
definition of a LIRG. However, given the head-on colli- 
sion between the Taffy galaxies, it is likely that they will 
eventually merge through dynamical friction. Thus the 
Taffy may move from its current pre-LIRG status to that 
of a LIRG or ULIRG if the two galaxies are drawn back 
together. Extended mid-IR H2 e mission is seen in ^ 70% 
of a sample of nearby ULIRGs (Higd on et al.l[2006l ). and 
shocks driven into th eir outer disks h ave been hypothe- 
sized but not proven ()Zakamskall201Clf ) . The existence of 
extended shocked H2 in the Taffy system 20 Myr after 
the initial collision suggests that H2 emission can persist 
for a time comparable with the dynamical crossing time. 
Thus given the even more chaotic nature of the gas-rich 
major mergers believed to explain many local ULIRGs, 
and their expected rapid coalescence rates, it seems plau- 
sible that the H2 emission seen in ULIRGs may result 
from shocked processes broadly similar to that seen in 
the Taffy. 



4.3. Warm H2 gas fraction 

One of the most remarkable features of the Taffy sys- 
tem is the fact that the peak H I concentration is in the 
bridge between the two galaxies (iCondon et al.lll993D . In 
contrast, the CO (1-0) map of iGao et al.l (|2003[ ) shows 
that the molecular gas has its highest concentrations 
in the galaxies. However, the molecular gas is present 
throughout the bridge, so that the total H2 mass in the 
bridge appears to be comparable to that of the galaxies. 

The estimated total H2 mass depends on the conversion 
factor X = iVHj /Ico , which is unknown. To estimate the 
warm gas fraction in the bridge, defined as the ratio of the 
warm H2 mass derived from the low-res IRS observations 
to the H2 mass derived from the CO observations, we first 
use the Galactic value of X = 2x10^° cm~^/( Kkms~^). 
The warm H2 mass is estimated to be 8 x 10* Mq based 
on the S(0) and S(l) lines (see Section [3^ . This gives a 
a warm gas fraction in the bridge of 0.09. We note that 
in the case of Stephan's Quintet, recent measurements 
with the IRAM 30 m telescope (jGuillard et al.l I2012D 
suggest comparable values along the main intergalactic 
shock MH2,warmA^H2 = 0.10-0.33 usiug the same Galac- 
tic X factor. 

However, the value of the X-factor in the Taffy bridge 
could be lower than the Galactic v a.lue by a factor of a 
few. Based on ^'^CO measurements. [Braine et "all ()20Q3f ) 
suggested that the X factor could be at least a factor 
of 4 lower than the Galactic value in the bridge region, 
which would increase the warm gas fraction by a factor 
of fo ur. 

.Zhu et al.l ()2007[ ) fit large velocity gradient models 
to CO line ratios and determined that X = 2.6 x 
10^^ cm~^/(K km s~^) in the bridge, which implies 
AfH2 = 1-3 X 10^ Mq and Mu^^^^^^/Mu^ = 0.7. There- 
fore, the warm H2 gas fraction in the bridge is poorly 
constrained, likely being in the range 0.09 — 0.7. This 
is on th e upper end of the distribution for the SINGS 
galaxies (jRoussel et al.ll2007l ) . Together with the absence 
of obvious source of UV heating, the bracketed range of 
warm H2 gas fraction suggests an active heat source from 
shocks. 

5. SUMMARY AND FUTURE WORK 

We have reported on the detection of powerful mid- 
IR emission lines of H2 in the Taffy galaxies. Our main 
results are: 

1. The mid-IR spectrum of the Taffy bridge bears a 
striking similarity to the spectrum of the group- 
wide shock in Stephan's Quintet. In particular, 
the L(H2)/L(PAH8/im) ratio is unusually high 
in the bridge, exceeding by an order of magni- 
tude that_Jound in SINGS star forming galax- 
ies (jRoussel et al.l 2007), while the galaxies UGC 
12914/5 are comparable to star forming galaxies. 

2. The warm H2 is distributed throughout the bridge 
between the two galaxies. This is similar to the 
H2 detected previously using sub-mm observations 
of CO. It is also comparable to the distribution of 
20 cm radio continuum and dust emission. 

3. Single-temperature fits to the two H2 lines mea- 
sured throughout the bridge indicate a nearly uni- 
form temperature of 150-175 K. We also estimate 
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a mass density ~5 x 10^ Mq kpc^^. We measure 
a warm H2 mass of 4.2 x 10* Mq in the observed 
parts of the bridge, and predict a total mass in 
the complete bridge of 9(+2/-5) x 10* Mq if the 
rest of the bridge has the same properties as that 
observed. 

4. The H2 surface luminosity in the Taffy bridge is 
more than a factor of two larger than the main 
shock region in SQ, probably because of different 
collision geometries and initial conditions. Since 
the H2 line cooling time is many orders of magni- 
tude shorter than the dynamical timescale for the 
bridge, in situ heating is required to explain the 
observed warm temperature of the bridge gas. The 
H2 gas is probably heated by turbulence and shocks 
produced during the collision, and also by contin- 
uing cloud collisions in the bridge. 

We have determined that the H2 molecule provides an 
important cooling channel for the warm gas in the Taffy 



bridge. Other molecules may also be important. Future 
observations (granted in Open Time 1) with the Herschel 
Space Observatory will allow us to further explore the 
thermodynamic state of the shocked ISM bridge. 
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Fig. 1. — IRS SH (a), LH (b), SL (c), and LL (d) slit positions overlaid on IRAC 8.0 fiiii image. Each module was nodded between two 
positions at each target, with the target at the 1/3 and 2/3 positions in the slit. A single nod position for each module is shown in blue for 
clarity. In (b), we also label the position we call the "bridge/arm" in the text. For the SL and LL modules, both orders are shown. North 
is up and east to the left. 




Fig. 2. — Spectral extraction regions A-U overlaid on MIPS 24 /jm. 
regions corresponds to one of the LL slit positions, which are labeled. 



The extraction regions are 10'' 15 X 10'' 15. Each row of extraction 
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Fig. 3. — High resolution spectra from tlie IRS SH (black) and LH (blue) modules. The SH spectra were scaled up to match the LH 
continua. The mean scaling factors for the two nods, along with the measured line fluxes, are presented in Table [T] 
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Low Resolution Spectra 
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Fig. 4. — Low resolution spectra from IRS LLl (blue), LL2 (black), SLl (green, where available) modules extracted from bridge regions 
A-G (position 1; left column) and H-O (position 2; right column). The SLl data have been scaled up to the same aperture size as the LLl 
assuming a uniform distribution over the aperture. Photometric data points from IRAC and MIPS are shown as red diamonds. 
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Fig. 5. — Low resolution spectra from IRS low resolution extraction regions P-S (position 3; left column) and T-U (position 4; right 
column). Symbols are as in Figure |4] 

Line Cross-Sections 



c 
3 



o 






c 




5 10 15 20 25 30 35 5 10 15 20 25 30 35 5 10 15 20 25 30 35 

pixels 



Fig. 6. — Cross sections through the LLl and LL2 slits at the wavelengths of [S IIl]33.48/.im, 30 /im continuum, H2 S(0) 28.22 ^m, 
[S IIl]18.71/im, 17.4 fim continuum, and H2 S(l) 17.03 ^m. The two continuum profiles give a sense of the the galaxy widths at wavelengths 
near the lines of interest. The positions are indicated in Figure [2] In position 2, the slit is slightly south of UGC 12915. These plots, 
particularly of the H2 S(l) line, demonstrate that H2 emission is extended across the full width of the Taffy bridge. 
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H2 S(0) 28.22 u 
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Fig. 7. — (a) 20 cin radio image Cond on et al.l II1993I) . with logarithmic color scale and contours from 0.1-7 mjy beam~^. (b) 20 cm 
contours are also shown overlaid on IRAC 8.0 ^m image (b), sparse maps of H2 S{1) 17.03 fim line (c), [S III18.71]/im line (d), H2 S{0) 
28.22 fim line (e), and [Si Il]34.82 fj,m line (f). The sparse maps of the S(0) and S(l) lines show that the H2 emission is extended across 
the entire bridge, compared to the emission from [S III] and [Si II], which have comparable wavelengths but do not show emission outside 
the galaxies. The units on plots (b)-(f) are MJy sr~^. 
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Fig. 8. — Excitation diagrams with multi-temperature fits of low-res regions C (green) and J (red). Tfie fit parameters are shown in 
Table [3] 
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Fig. 10. — Total luminosity in H2 0-0 S(0) and S(l) lines relative to PAH emission in the IRAC 8.0 /^m band. Extraction regions which 
fall on the disks or their nuclei have filled symbols. The PAH emission was determined using rectangular apertures matched to the spectral 
extraction apertures in the IRAC 8.0 fiicn band, which was Gaussian-convolved to match the resolution of the MIPS 24 fira band and 
corrected for stellar contamination (see text). Luminosities at 24 /im were determined using MIPS 24 fim photometry, except for the high 
resolution Taffy apertures, for whichjm average of the spectrum o ver 22.5-25.0 /jm was used (see text). For comp arison, we sh ow galaxies 
from the SINGS sample I Rouss ellrt a L. 200l!) . Arp 143 knot G [Be irao et al.ll200g) , and the SQ shock sub-region llCluver et al. 2010). 
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Fast Shock MAPPINGS models - B = 5/iG 
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Fig. 12.— MAPPINGS shock velocity models HAllen et al.|[200a ) for magnetic field B = S/xG. The dotted line indicates [Ne IIl]/[Ne II] 
= 0.62, as measured in the bridge using the IRS SH module. The curves correspond to different values for the pre-shock density. 
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TABLE 
H2 Line Fluxes (IC 


1 

-1^ W M-2) 










Region 


RA 


Dec 


H2 0-0 S(0) 


Ha 0-0 S(l) 


Ha 0-0 S(2) 


Ha 0-0 S(3) 


Scaling'^ 




(J2000.0) 


(J2000.0) 


A28.22^m 


A17.03^m 


A12.28Atm 


X9.66fim 




High Resolution'' 


UGC 12915 nucleus 


01 41.89 


+23 29 44.0 


3.17 + 0.25 


10.41 ± 0.19 


3.02 ± 0.09 




SH; 1.42 


UGC 12914 nucleus 


01 38.09 


+23 29 03.3 


1.57 ± 0.21 


6.57 ± 0.14 


2.30 ± 0.10 




SH; 1.35 


UGC 12914 S 


01 39.02 


+23 28 43.0 


2.19 ± 0.27 


4.89 ± 0.10 


1.51 ± 0.08 




SH; 1.42 


UGC 12914 N 


01 37.37 


+23 29 18.5 


1.06 ± 0.13 


2.53 ± 0.12 


0.72 ± 0.10 




SH; 1.93 


bridge/arm 


01 38.89 


+23 29 30.0 


1.26 ± 0.13 


3.61 ± 0.09 


1.25 ± 0.10 




SH; 2.81 


Low Resolution'^ 


A 


01 37.62 


+23 29 18.7 


0.50 ± 0.16 


2.10 ± 0.24 








B 


01 38.29 


+23 29 22.9 


0.71 ± 0.08 


3.21 ± 0.27 










C 


01 38.97 


+23 29 27.0 


0.63 ± 0.04 


2.65 ± 0.16 










C (PAHFIT)'* 


01 38.97 


+23 29 27.0 


0.61 ± 0.09 


2.74 ± 0.20 


0.40+0.03 


0.61+0.04 


SLl, 


SL2; 3 


D 


01 39.64 


+23 29 31.1 


0.68 ± 0.08 


3.18 ± 0.17 










E 


01 40.32 


+23 29 35.2 


0.95 ± 0.12 


4.67 ± 0.16 










F 


01 40.99 


+23 29 39.3 


1.33 ± 0.03 


6.88 ± 0.18 










G 


01 42.00 


+23 29 45.5 


1.03 ± 0.28 


6.31 ± 0.22 










H 


01 38.06 


+23 29 04.8 


0.61 ± 0.16 


4.94 ± 0.25 










I 


01 38.74 


+23 29 08.9 


0.46 ± 0.04 


3.18 ± 0.12 










J 


01 39.41 


+23 29 13.0 


0.51 ± 0.09 


2.72 ± 0.13 










J (PAHFIT)d 


01 39.41 


+23 29 13.0 


0.72 ± 0.06 


2.90 ± 0.16 


0.37+0.02 


0.56+0.03 


SLl, 


SL2; 3 


K 


01 40.09 


+23 29 17.1 


0.72 ± 0.15 


2.83 ± 0.21 










L 


01 40.77 


+23 29 21.3 


0.48 ± 0.08 


2.78 ± 0.21 










M 


01 41.44 


+23 29 25.4 


0.30 ± 0.06 


1.41 ± 0.12 










N 


01 42.12 


+23 29 29.5 


< 0.41= 


1.02 ± 0.22 










O 


01 42.79 


+23 29 33.6 


< 0.25° 


0.65 ± 0.12 










P 


01 39.15 


+23 28 43.6 


0.65 ± 0.18 


2.62 ± 0.21 










Q 


01 40.16 


+23 28 49.8 


0.50 ± 0.13 


0.86 ± 0.16 










R 


01 40.84 


+23 28 53.9 


0.36 ± 0.10 


0.69 ± 0.25 










S 


01 41.51 


+23 28 58.0 


< 0.16= 


0.93 ± 0.12 










T 


01 39.60 


+23 28 29.7 


< 0.25= 


< 0.46= 










U 


01 40.61 


+23 28 35.9 


< 0.27= 


< 0.46= 











'^ Due the different angular size of the IRS modules, some of the spectra were rescaled so that the continua matched. The 

module and scaling factor are indicated here. 

^ Ha S(0) measured with Spitzer IRS LH module; S(l) and S(2) lines measured with IRS SH module. All fluxes were 

measured using SMART. 

= Ha S(0) line measured in LLl, S{1) in LL2, S(2) and S(3) where available in SLl. Low resolution bridge region apertures 

were squares 10" 15 X 10" 15. Fluxes were measured using SMART except where otherwise noted. 

'^ Line fluxes measured by fltting the full spectrum using PAHFIT, which requires a smooth continuum. In region J, SLl 

and SL2 were both scaled up by a factor of 3 due to their smaller area. Then SLl, SL2, and LL2 were scaled up by a factor 

1.86 to match the LLl continuum. Following the fltting with PAHFIT, the SLl, SL2, and LL2 were scaled down by 1.86. 

No rescaling was required in region C. 

= 3(T upper limit estimated from the RMS and expected line width. 



TABLE 2 
Fine structure line fluxes for high resolution data (lO"'-^ W m" 



Region [Ne II] [Ne III] [S III] 

A12.8lMm A15.56Afm A18.7lMm 



[FeIl] + [OlV] 
A25.99^m+A25.89Mm 



[S III] 
A33.48/im 



[Si II] 
A34.82/im 



UGC 12915 nuc 
UGC 12914 nuc 
UGC 12914 S 
UGC 12914 N 
bridge/arm 



27.94+1.08 
4.26+0.28 
9.14+0.39 
6.35+0.32 
0.58+0.06 



3.79+0.18 
1.52+0.11 
1.42+0.07 
1.53+0.07 
0.36+0.14 



10.25+0.45 
1.79+0.15 
5.96+0.17 
4.20+0.20 
0.68+0.07 



1.91+0.53 
0.65+0.15 
0.33+0.10 
0.38+0.14 
< 0.23=' 



18.72+1.61 
2.89+0.59 
9.95+0.73 
7.93+0.47 
0.58+0.23 



28.06+0.68 
8.34+0.37 
10.78+0.43 
9.41+0.54 
1.63+0.28 



Note. — [Ne II], [Ne III], and [S IIl]18.71/xm measured with IRS SH module; [Fe Il] + [0 IV], [S IIl]33.48/im, 
and [Si II] measured with IRS LH module. The SH line fluxes have been scaled up using the same factors as in 
Table [1] 
^ 3a upper limit estimated from the RMS and expected line width. 
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TABLE 3 
Ht Properties of Regions A-R 



Region 



Temperature'' Equilibrium o/p° ^11-2 Sjja 

(K) (lO^Omolceules cm^^) (iqS Mq kpe^^) 



A 


160 


±10) 


2.6 (±0.1) 


2.7 (+1.0/-0.9) 


4.4 (+1.4/-1.2) 


B 


163 


±10) 


2.6 (±0.1) 


3.7 (+1.1/-0.9) 


5.9 (+1.5/-1.2) 


C 


160 


±5) 


2.6 (±0.1) 


3.4 (±0.4) 


5.5 (+0.8/-0.7) 


C-multi Ti 


160 


±5) 


2.6 (±0.1) 


3.2 (±0.4) 


5.0 (+0.8/-0.7) 


C-multi T2 


488 


+10/-70) 


3.0 


0.010 (+0.01/-0.001) 


0.020 (+0.020/-0.001) 


D 


164 


+8/-7) 


2.6 (±0.1) 


3.5 (+0.8/-0.7) 


5.6 (+1.4/-1.2) 


E 


167 


±2) 


2.7 (±0.1) 


4.6 (±0.2) 


7.5 (±0.3) 


F 


169 


±2) 


2.5 (±0.1) 


6.4 (+0.4/-0.3) 


10.2 (+0.6/-0.4) 


G 


178 


±13) 


2.7 (±0.1) 


4.5 (+2.2/-1.6) 


7.2 (+2.8/-2.5) 


H 


195 


±12) 


2.9 (±0.1) 


2.2 (+1.0/-0.8) 


3.6 (+1.7/-0.7) 


I 


185 


±5) 


2.8 (±0.1) 


1.9 (+0.5/-0.4) 


3.0 (+0.5/-0.4) 


J 


165 


±3) 


2.6 (±0.1) 


2.6 (+0.4/-0.1) 


4.2 (+2.0/-0.2) 


J-multi Ti 


155 


±3) 


2.5 (±0.1) 


4.1 (+0.4/-0.1) 


6.6 (+2.0/-0.2) 


J-multi T2 


433 


±3) 


3.0 


0.020 (±0.001) 


0.020 (+0.001/-0.001) 


K 


157 


±10) 


2.6 (±0.1) 


4.1 (±1.8) 


6.6 (+2.7/-0.1) 


L 


175 


±11) 


2.7 (±0.2) 


2.1 (±0.7) 


3.4 (+1.2/-1.2) 


M 

N 


165 


±11) 


2.6 (±0.1) 


1.5 (±0.6) 


2.4 (+1.1/-0.8) 



P 


158 


±13) 


2.6 (±0.2) 


3.6 (±2.0) 


5.8 (+3.2/-2.5) 


Q 


130 


+14/-11) 


2.2 (±0.2) 


4.6 (+3.0/-2.0) 


7.3 (+5.0/-3.0) 


R 


133 


+10/-12) 


2.3 (+0.1/-0.3) 


3.1 (+1.0/-0.2) 


5.0 (+1.8/-1.0) 


UGC 12915 nuc Ti 


143 


+5/-6) 


2.4 (±0.1) 


22 (+4/-3 ) 


36 (+6/-5) 


UGC 12915 nuc T2 


696 


+30/-23) 


3.0 


0.07 (±0.01) 


0.11 (±0.02) 


UGC 12914 nuc Ti 


151 


+8/-6) 


2.5 (±0.1) 


9.6 (+3.0/-2.0) 


15 (+4/-3) 


UGC 12914 nuc T2 


793 


+10/-6) 


3.0 


0.040 (±0.005) 


0.060 (±0.008) 


UGC 12914 N Ti 


135 


+8/-7) 


2.3 (±0.1) 


8.6 (+3/-2) 


14 (+4/-3) 


UGC 12914 N T2 


1500 


(±10) 


3.0 


0.010 (±0.005) 


0.020 (±0.008) 


UGC 12914 S Ti 


132 


+4/-3) 


2.8 (±0.1) 


19 (+4/-3) 


30 (+8/- 7) 


UGC 12914 S T2 


989 


±40) 


3.0 


0.020 (±0.005 ) 


0.040 (±0.008) 


Hi res bridge/arm single-temp fit 


157 


+23/-14) 


2.6 (±0.2) 


2.3 (+1.7/-1.3) 


3.6 (+3.0/-1.9) 


Hi res bridge/arm 2-temp fit Ti 


103 


+27/-3) 


1.7 (+0.5/-0.2) 


10 (+3/-4) 


15 (+5/-6) 


Hi res bridge/arm 2-temp fit T2 


310 


±30) 


3.0 


0.12 (-I-0.5/-0.7) 


0.2 (+0.05/-0.12) 



'^ Uncertainties in all derived properties are formal uncertainties in fitting and do not include many possible systematic 

effects 

^ Equilibrium ortho-to-para ratio for H2 molecules is assumed. Note that this can be significantly less than 3 for T < 300 K, 

though at high temperatures o/p = 3 with no formal uncertainty apart from that in the temperature determination. For 

Deviations from thermal equilibrium might be possible in shocks and this would further add to the uncertainties in the 

derived properties. 



